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Abstract 



The largely open problem of scaling laws in fully developed turbulence is discussed, 
with the stress put on similarities and differences with scaling in field theory. A soluble 
£> ■ model of the passive advection is examined in more detail in order to illustrate the 

principal ideas. 
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1 LANGEVIN VERSUS NAVIER- STOKES 



Many dynamical problems in physics may be described by evolution equations of the type 



X ■ dt$= - F($) + / (1.1) 

h ; 

where 3>(£, x) represents local densities of physical quantities, i ? ( < 5) is their nonlinear func- 
tional and / stands for an external source. We shall be interested in the situations where 
the source / is random. For concreteness, we shall assume it Gaussian with mean zero and 
covariance 

</(f,x)/(*,y)> = 5(t-s)C(^) (1.2) 

with L determining the scale on which f(t, x) are correlated. An example of such an 
evolution law is provided by the Langevin equation describing the approach to equilibrium 
in systems of statistical mechanics or field theory Q . In this case the nonlinearity is of the 
gradient type: 
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with S($) a local functional, e.g. = \ J(V$) 2 + \m 2 J $ 2 + A / <£> 4 in the $ 4 field 

theory, and with L small so that C(x/L) is close to the delta-function <5(x) and regulates 
the theory on short distances < L. 

Augmented by an initial condition <I>(io,x) (and, eventually, boundary conditions), the 



solution of Eq. (LI) should define random field <l?(i,x). We are interested in the behavior 



of its correlation functions given by the mean values 

i 

Among the basic questions one may ask are the following ones: 

1. Do the correlation functions become stationary (i.e. dependent only on time 
differences) when to — > — oo? If so, are the stationary correlators unique (independent of 
the initial condition)? 

2. Do they obey scaling laws? 

For the field theory case these questions are well studied with the use of powerful analytic 
tools as perturbative expansions and renormalization group and by numerical analysis 
(Monte Carlo simulations). The stationary correlators describe possibly different phases 
of the system. Universal (i.e. independent of the cutoff C) scaling laws of the type 

<Q(*(i,x))Q($(i,y))) ~ |x-y[-fc (1.5) 

for |x — y 3> L with Q some local functions of f> emerge at the points of the 2 nd order 
phase transitions. 

On the opposite pole of the field theoretic the hydrodynamical examples of the 

evolution equation fll.l|) . The best known of those is the Navier-Stokes equation 

d t v = - P(v • V)v + i/Av + f (1.6) 

for the incompressible (V • v = 0) velocity field v(t, x) , with P standing for the orthogonal 
projection on such vector fields, v denotes the viscosity and f is the external force which 
induces the fluid motion. In the fully developed turbulence one is interested in the regime 
where the stirring forces act on large distances (like the convective forces on scales of 
kilometers in the atmosphere) and we observe quite complicated (turbulent) motions on 
shorter distances down to scales on which the dissipative term cx v becomes important ( 
~ milimeters in the atmosphere). It is believed that the large scale details should not be 
essential for the statistics of the flow in this intermediate regime called the "inertial range" . 
It is therefore common to model the stirring forces by a random Gaussian process with 
mean zero and covariance 

(r(i,x)/( S ,y)> = 6{t-8)C#(?jZ) (1.7) 

with d a C al3 = 0. L denotes now the large "integral scale" on which the random forces act. 
Note that, unlike for field theory, in this case the covariance C(x/L) is close to a constant, 
i.e. to a delta-function in the wavenumber space and not in the position space. Such regime 
in field theory would correspond to distances shorter than the ultraviolet cutoff with the 
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behavior strongly dependent on the detailed form of the regularization. Another (related) 



difference is that in Eq. (1.6) the nonlinear term is not of a gradient type. Finally, the 
projection P renders it non-local which is another complication. All these differences make 
the Navier-Stokes problem ( |1.6[) quite different from that posed by the Langevin equation 



and resistant to the methods employed successfully in the study of the latter. 

2 KOLMOGOROV THEORY 

The first major attempt to obtain universal scaling laws for the inertial range correlators 
is due to Kolmogorov Q. Assuming the existence of homogeneous (i.e. translationally 
invariant) stationary correlators of velocities, one deduces the following relation at equal 
times 

I V x • ( (v(x) - v(y)) 2 (v(x) - v(y)) ) + \ tr C(^) = u{ Vv(x) • Vv(y) ) . (2.1) 



Eq. ( |2.lD is obtained the following way. First, using Eq. fll.q ), we write 

t+st 

v(t + (5t,x) = v(i,x) + [-P(v- V)v(i,x) + uAv(t,x)]5t + J f{s)ds + 0(5t 2 ) (2.2) 

V 

and we instert this expansion into the equal-time 2-point function obtaining 

(v(t + 6t,XL)v(t + 6t,y)) = (v(t,x)v(t,y)) 

-(P(v- V)v(t,x).v(t,y))ft- (v(t,x)-P(vV)v(t,y)<ft 

+ v{Av(t,x) -v(t,y))5t + v(v(t,x) ■ Av(t,y))5t + tr C(^-)St + 0(5t 2 ) . (2.3) 

Vanishing of the 0(5t) terms produces Eq. (|2.1[) . 

Under the limit y — > x for positive v, Eq. ( |2.1| ) becomes the identity 

itrC(O) = (KVv) 2 ) = e (2.4) 

which expresses the energy balance: in the stationary state, the mean energy injection rate 
is equal to the mean rate of energy dissipation e. On the other hand, performing the limit 
v — > for x^y one obtains 

-iv x -((v(x)-v(y)) 2 (v(x)-v(y))} = ^trC(^). (2.5) 

For |x — y| <C L the right hand side is approximately constant and equal to ^-tr C(0), i.e. 
to e. Assuming also isotropy (rotational invariance) of the stationary state one may then 
infer the form of the 3-point function in the inertial range: 

(K(x)-^(0))(^(x)-^(0))(^(x)-^(0))) - --^^(^/V + cycL). (2.6) 

The latter implies for the 3-point function of the component of (v(x) — v(0)) paralel to x 
the relation 

<M*>-*°»?> a -55TT2) £ " |x| <2 ' 7) 
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known as the Kolmogorov "| law" (| is the value of the coeffitient on the right hand side 
in 3 dimensions). 

Under natural assumptions about the stationary state one may deduce a stronger ver- 
sion (sometimes called refined similarity) of the above relation which takes the form of the 
operator product expansion for the v — > limit of the dissipation operator e = y(Vv) 2 : 



e(x) = -j hm V x ■ [(v(x) - v(y)) 2 (v(x) - v(y))] 



(2.8) 

holding inside the expectations in the v — > limit. Relation fl2.8|) expresses the dissipative 
anomaly: the dissipation e whose definition involves a factor of v does not vanish when 
v -> 0. 

Kolmogorov postulated § that the scaling of general velocity correlators in the inertial 
range should be determined by universal relations involving only the distances and the 
mean dissipation rate e. Such postulate leads to the scaling laws 

((v(x) - v(0))™) oc e n/3 |x| n/3 (2.9) 

for the n-point "structure functions" of velocity generalizing the (essentially rigorous) result 
fl2.7| ) about the 3-point function. The right hand side of relation ( |2.9| ) is the only expression 
built from e and |x| with the right dimension. 

The power law fits oc |x|^ n for the structure functions measured in experiments and 
in numerical simulations lead to the values of the exponents slightly different from the 
Kolmogorov prediction ( n = n/3 for n ^ 3. One obtains |§ (2 — -70, £4 = 1.28, ^6 — 
1.77, ^8 — 2.23. The discrepencies indicate that the random variables (v(x) — v(0))., 
are non-Gaussian for small x with the probability distribution functions decaying slower 
than in the normal distribution. Such a slow decay signals the phenomenon of frequent 
occurence of large deviations from the mean values called " intermittency" . There exist 
many phenomenological models of intermittency of the inertial range velocity differences 
based on the idea that the turbulent activity is carried by a fraction of degrees of freedom 
with a self-similar (" multi- fractal" ) structure An explanation of the mechanism behind 
the observed intermittency starting from the first principles (i.e. from the Navier-Stokes 
equation) is, however, still missing and constitutes the main open fundamental problem of 
the fully developed turbulence. 



3 KRAICHNAN MODEL OF PASSIVE ADVECTION 

Recently some progress has been achieved fH,||,f7j in understanding the origin of inter- 
mittency in a simple model describing advection of a scalar quantity (temperature 
T(i,x)) by a random velocity field v(i,x). The evolution of the temperature is described 
by the equation 

d t T = - (v • V)T + kAT + f (3.1) 

where k denotes the molecular diffusivity and / is the external source which we shall take 
random Gaussian with mean zero and covariance (|h^). Following Kraichnan H, we shall 
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assume that v(t,x) is also a Gaussian process, independent of /, with mean zero and 
covariance 

(« a (t,x)^(s,y)> = 6(t-8)(D o S°0-d«f ) (x-y)) (3.2) 

with Dq a constant, cP^(x) oc |x|^ for small |x| and with d a d a ^ = in order to assure 
the incompressibility. Note the scaling of the 2n-point function of velocity differences with 
power n£ of the distance. The Komogorov scaling corresponds to £ = | (the temporal delta- 
function appears to have dimension length^ -1 ). The time decorrelation of the velocities 
is not, however, a very realistic assumption. In the Kraichnan model £ is treated as a 
parameter running from to 2. 

Writing 

t+St 

T(t + 5t,x) = T(i,x)- J [(v-V)T(s,x) - /(s, x)] ds + nAT(t, x)6t + 0(5t 2 ) , (3.3) 



we obtain the analogue of the relation (2.1) for the stationary state of the scalar (the latter 
may be shown to exist and to be independent of the initial condition decaying at spatial 
infinity) : 

- i^(x-y)^cy (T(x)T(y)) + |C(^) = k ( VT(x) ■ VT(y) ) . (3.4) 
Letting in Eq. ( |3.4|) y — > x for k > produces, as before, the energy balance: 

JC(0) = (k(VT) 2 } = 6. (3.5) 
On the other hand, taking k — > for x/y results in the equation analogous to (|2.5|): 



i^(x-y)^^(T(x)T(y)) = Jc(^) (3.6) 
which may be easily solved exactly for the 2-point function of the scalar giving 

(T(x)T(O)} = A c L 2 ~t - const. e|x| 2 -« + 0(L~ 2 ) (3.7) 

or 

( (T(x) - T(0)) 2 ) ^ const, e |x| 2 ~^ (3.8) 

for |x| <C L. This is an analogue of the Kolmogorov | law. It may be strengthen to the 
operator product expansion for the dissipation operator |1C|],|15] e = k(VT) 2 



e(x) = i lim Jpf > {x-y)d a T(x)d fl T(y) 



2 y^ x 



(3.9) 

K=0 



valid inside the expectations in the limit k — > 0. Eq. (3.9) expresses the dissipative anomaly 
in the Kraichnan model, analogous to the dissipative anomaly fl2,8| ) for the Navier-Stokes 
case. 
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The natural question arises whether the higher structure functions of the scalar ( (T(x)— 
T(0)) 2n ) = 52n(x) scale with powers n(2 — £) as the dimensional analysis would suggest 
(Corrsin's analogue [pi of the Kolmogorov theory). The answer is no. Experiments show 



that the scalar differences display higher intermittency than that of the velocities [12|. 
Although, by assumption, in the Kraichnan model there is no intermittency in the distri- 
bution of the velocity differences, numerical studies Jl3|],|l4j indicate strong intermittency 
of the scalar differences signaled by anomalous (i.e. 7^ n(2 — £)) scaling exponents. Unlike 
in the Navier-Stokes case, we have now some analytic understanding of this phenomenon, 
although still incomplete and controvertial ||,[|l4j|,[16|. 



The simplifying feature of the Kraichnan model is that the insertion of expansion ([: 
into the higher point functions J r 2n(x) = {T(t, xi, . . . , T(t, X2 n )) leads to a system of (Hopf) 
equations which close: 

-M2n^2n(x) = £ ^2n-2^l, ■ • » , X 2n ) C(^) (3.10) 

i<j i j 

where Ai n are differential operators 

n 

M n = d af3 ( Xi - Xj ) d x? d x p - kJ2 A Xi ■ (3.11) 

l<i<j<n i=l 

In principle, the above equations permit to determine uniquely the stationary higher-point 
correlators of the scalar iteratively by inverting the positive elliptic operators M.2n- By 
analyzing these operators whose symbols loose strict positivity when k goes to zero, it was 
argued in |5[,p5[ that at least for small £ 



1. lim J-2n(x) exists and is finite, 

2. ^2n(x) = A c>2n L^ ^ 2 °„(x) + O(L- 2 +°(0) + . . . a t « = 

where p2n = 2 "i" ^ + C(£ 2 )> is a scaling zero mode of the n = version M.2 n of 
the operator M 2n - M° 2n J^ n = 0, ^ n (Ax) = \^~0-p2n jrOjx) . The ...-terms do not 
depend on at least one of the vectors x» and do not contribute to the correlators of scalar 
differences. 

3. 5 2n (x) oc L P2n | x [ n ( 2 -5)-P2n at k = and for |x| < L. 

The last relation, a simple consequence of the second one, shows appearence of intermittent 
exponents at least for small £ (p2n is their anomalous part and it is positive starting from 
the 4-point function). A similar analysis, consistent with the above one, has been done for 
large space dimensions |6||,[p7|. 



The above results about the "zero-mode dominance" of the correlators of the scalar dif- 
ferences show what degree of universality one may expect in the scaling laws of intermittent 
quantities: the amplitudes Ac. 2n in front of the dominant term depend on the shape of 
the covariance C i.e. on the details of the large scale stirring. But the zero modes of the 
dominant terms (and their scaling exponents) do not. In field theory, small-scale univer- 
sality of the critical behavior finds its explication in the renormalization group analysis. 
Similarly, in the Kraichnan model there exists a renormalization group explanation of the 
observed long scale universality [18]. The renormalization group transformations eliminate 
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subsequently the long scale degrees of freedom. In a sense, they consist of looking at the 
system by stronger and stronger magnifying glass so that the long distance details are lost 
from sight. The eliminated degrees of freedom induce an effective source for the remaining 
ones |L8|| . Whereas such an "inverse renormalization group" analysis may be implemented 
for more complicated turbulent systems rests an open problem. 



4 DYNAMICS OF LAGRANGIAN TRAJECTORIES 

What is the source of the zero mode dominance of the inertial range correlators of the 



scalar differences? In absence of the diffusion term in Eq. (3.1) the scalar density is given 
by the integral 

t 

T(i,x) = / f(s;y(s;t,x))ds (4.1) 



where y(s; t, x) describes the Lagrangian trajectory, i.e. the solution of the equation 

^=v( S ,y), (4.2) 
as 

passing at time t through point x (for concretness, we have assumed the vanishing initial 
condition for T at t = — oo). The Lagrangian trajectories describe the flow of the fluid 
elements. If the velocities are random, so are Lagrangian trajectories and we may ask the 
question about their joint probability distributions. Let P n (t, x; s, y) denote the probability 
that n Lagrangian trajectories starting at time s at points (yi, . . . ,y n ) = y pass at time 
t > s through points (xi,...,x n ) = x. These probabilities may be computed for the 
time-decorrelated Gaussian velocities. They appear to be given by the heat kernels of the 
singular elliptic operators M.^ : 

Pn(t,x;s,y) = e-(*- s ^(x;y) (4.3) 

(more exactly, this is true after averaging over the simultaneous translations of the initial 
points, i.e. for the probabilities of relative positions of the trajectories). 

From the form of the operators A4^ we see that the (relative positions of) n La- 
grangian trajectories undergo a diffusion process with distance-dependent diffusion co- 
efficients. When two trajectories are close, the corresponding diffusion coefficient vanishes 
as the distance to power £ slowing down the diffusive separation of the trajectories. When 
the trajectories eventually separate, the diffusion coefficient grows speeding up further 
separation. The result is a super diffusive large time asymptotics: 

/ /(x) P n {t,x;0,y)dx oc t& (4.4) 

for a generic (translationally invariant) scaling function /(x) of scaling dimension a > 
(e.g. for J2i j |xj — Xj| 2 with the scaling dimension 2). Note the faster than diffusive growth 
in time for £ > 0. There are, however, exceptions from this generic behavior. In particular, 
if / is a scaling zero mode of M.^ then 

J /(x) P„(t,x; 0,y) dx = J /(x) e- tM °(x;y) dx = const. = /(y) . (4.5) 
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It can be shown [18] that each zero mode fo of A4^ of scaling dimension <to > generates 
descendent slow collective modes f p , p=l,2,. . . , of scaling dimensions a p = <to +p{2 — £) 
for which 

J / P (x)Pn(*,x;0,y)dx oc t P, (4.6) 



i.e. grows slower (if do > 0) than in ( |4.4J ) . The descendants satisfy the chain of equations 
M-nfp = fp-i- The structure with towers of descendants over the primary zero modes 
resembles that in systems with infinite symmetries and may suggest presence of hidden 
symmetries in the Kraichnan model. 

The slow modes f p appear in the asymptotic expansion 

oo 

P n (t,x;0 5 y/L) = ]T ^iT^^i) f p (y) (4.7) 

zero modes p=Q 

fo 

valid for large L and describing the behavior of the trajectories starting close to each other. 
Since P n (i,x; 0, y/L) = P n (t, y/L; 0, x) , expansion flL7[ ) describes also probabilities that 
the trajectories approach each other after time t. By simple rescalings and the use of 
expansion (p~7|), one obtains 

/ /(x) P„(i,x; 0,y) dx = t^l J /(x) P„(l,x; 0,y/t^) dx 

= E fp(y) //(x)s p (i,x)dx. (4-8) 

zero modes 53— Q 

/o 

For generic /, the dominant term comes from the constant zero mode and is proportional 
to t 2 ~z . However for / equal to one of the slow modes, the leading contribution is given by 
the further terms in the expansion due to the orthogonality relations between / p 's and g p J s. 



The behaviors (|4.4|j4.5| , |4.6D result. The distributions P n (t,x; 0,y) enter the expressions for 
scalar correlators. For those of scalar differences, the contributions to expansion ( [4.7] ) of 
modes which do not depend on all y^ (like the constant mode) drop out leading to the 
dominance by non-trivial zero modes. 



The asymptotic expansion ( |4.7| ), governed by the slow modes f p , displays another im- 
portant feature of the Lagrangian trajectories. Since P n (t,x; 0,y) is a joint probability 
distribution of the endpoints of n Lag rangian trajectories starting at points yi, . . . ,yn, we 
should expect that it becomes concentrated on the diagonal xi = • • • = x n when the initial 
points yi, . . . ,y n tend to each other. But this is not the case as lim P n (t,x; 0,y/L) = 

L^oo — 

Pn(t, x; 0, 0) is a finite function of x. It is the constant zero mode contribution to expansion 
Q. The whole expansion describes how exactly P n (t,x; 0,y/L) fails to concentrate on 
the diagonal when L — > 00. Similarly, P n (i, 0; 0, y) is a non-singular function of y showing 
that the probability that the trajectories collapse after time t to a single point is finite, 
contradicting the uniqueness of the Lagrangian trajectories passing through a given point. 
The solution of the paradox is as follows. The typical velocities in the ensemble that we 
consider have rough spatial behavior (and even rougher time behavior). As functions of 
x they are (essentially) Holder continuous with exponent j. But for such velocities, the 
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equation for the Lagrangian trajectories (fL^) does not have a unique solution, given the 
initial position. As a result, the Lagrangian trajectories loose deterministic character for 
a fixed velocity realization. Nevertheless, one may still talk about probability distribution 
fn(i)X) 0,y|v) of their final points whose average over v gives P n (t,x; 0, y). In a more 
realistic description which takes into account a smoothing of the typical velocities at very 
short viscous scale, the same effect is due to the sensitive dependence of the now determin- 
istic Lagrangian trajectories on the initial conditions, within the viscous scale, signaled by 
the positivity of the Lyapunov exponent [jigf ]. 

Summarizing, intermittency in the Kraichnan model of the passive advection appears 
to be due to the slow collective modes in the otherwise super diffusive stochastic Lagrangian 
flow. The presence of such modes is closely related to the breakdown of the deterministic 
character of Lagrangian trajectories for the fixed velocity configuration at high Reynolds 
numbers, due to the sensitive dependence of the trajectories on the initial conditions within 
the viscous scale. We expect both phenomena to be present also in more realistic turbulent 
velocity distributions and to be still responsible for the anomalous scaling and intermit- 
tency. 
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